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C/3 ' Surface waves can propagate on the planar interface of a linear electro-optic (EO) material 

O ' and an isotropic dielectric material, for restricted ranges of the orientation angles of the EO 

lyj . material and the refractive index of the isotropic material. These ranges can be controlled by 

^-^' the application of a dc electric field, and depend on both the magnitude and the direction of 

rS^ [ the dc field. Thus, surface-wave propagation can be electrically controlled by exploiting the 

CiH. Pockets effect. 
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^^ '. 1 Introduction 
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Wave propagation localized to a planar metal-dielectric interface has a long history (Zenneck, 
1907; Agranovich and Mills, 1982; Boardman, 1982; Homola et al., 1999; Matveeva et al., 2005). 
f~^ , In contrast, wave propagation localized to the planar interface of an isotropic dielectric material 

and a uniaxial dielectric material was found possible by D'yakonov only in 1988 (D'yakonov, 1988; 
Averkiev and Dyakonov, 1990). Since then, researchers have explored surface-wave propagation 
(SWP) on increasingly complex systems of bimaterial interfaces such as biaxial-isotropic, uniaxial- 
uniaxial, and biaxial-biaxial (Walker et al., 1998; Darinskii 2001, Wong et al., 2005; Polo et al., 
2006, 2007a, b; Nelatury et al., 2007). In all such investigations, essentially, one looks for the selected 
angular regimes of the propagation direction wherein certain dispersion conditions are met. For 
complex systems, the angular regimes are very narrow and depend highly on the crystallographic 
symmetries of the two materials. 

Our motivation for this paper is to show how one might exercise control on the angular regimes 
of SWP. Whereas temperature and pressure may be altered to control SWP, the application of 
an external field is expected to provide dynamic control. In particular, the electro-optic effect 
refers to changes in optical properties by the application of a low-frequency or dc electric field. For 



instance, an optically isotropic crystal (possessing either the 43?Ti or 23 point group symmetry) 
upon exposure to a dc electric field turns birefringent (Cook, 1996; Lakhtakia, 2006a). 

The Pockels effect is a linear electro-optic (EO) effect, whereby the modification of the in- 
verse (optical) permittivity matrix by the dc electric field is quantified through 18 electro-optic 
coefficients, not all of which may be independent of each other, depending on the point group 
symmetry (Boyd, 1992). The opportunities offered by the Pockels effect for tuning the optical 
response characteristics of materials have recently been highlighted for photonic band-gap engi- 
neering (Lakhtakia, 2006b; Li et al., 2007) and composite materials (Lakhtakia and Mackay, 2007; 
Mackay and Lakhtakia, 2007). These publications suggest that, although the changes in the optical 
permittivity matrix are typically small, the effect on SWP could be significant due to the extreme 
sensitivity of surface waves on the dielectric constitutive properties. Let us note here an earlier 
study wherein the quadratic EO effect named after Kerr was shown to offer control over SWP 
(Torner et al., 1993). 

This paper introduces the influence of the Pockels effect on SWP. Although the theoretical 
treatment is general, numerical results are presented only for a specific EO material: potassium 
niobate (Zgonik et al., 1993). The remainder of the paper is organized as follows: Section 2 
provides a description of a canonical boundary-value problem, the optical permittivity matrix of 
an EO crystal, and the derivation of the dispersion relations for SWP. In Section 3 numerical 
results are furnished, and our conclusions are distilled in Section 4. A note on notation: vectors 
are underlined, matrixes are decorated with an overbar, and the Cartesian unit vectors are denoted 
as u^, Uy, and u^. All field quantities are assumed to have an exp(— iwi) time-dependence. 

2 Theory 

Let the plane of SWP be a bimaterial interface. The half-space z < is filled with a homogeneous, 
isotropic, dielectric material with an optical refractive index denoted by Ug. The half-space z > 
is filled with a homogeneous, linear, EO material, whose optical relative permittivity matrix is 
stated as (Lakhtakia and Reyes 2006a, b) 

erel = Sz (V') • Ryix) ' ^PE ' Ryix) ' Sz ("V') ■ (1) 

Incorporating the Pockels effect due to an arbitrarily oriented but uniform dc electric field E_ '', 
the matrix epE is given by 

/ M(. ^(0)^-3 pdc\ MMx^S ^ pdc ^(0)^(0)^-3 ^ pdc \ 

PDE. ~ ^(0)^(0) Y-3 ^ 77.dc ^(0)/-, ^(0)y-3 ^ T^dc\ MMs^-i ^ wdc 

(2) 



correct to the first order in \E '^l, where 



El' 1 ^Ry{x)-Sz{-i^)- I i?f I , (3) 

\e: 



dc 




^12 3 ^^^ th^ principal relative permittivity scalars in the optical regime, whereas rjK (with 1 < 
J < 6 and I < K < 3) are the EO coefficients. The EO material can be isotropic, uniaxial, 
or biaxial, depending on the relative values oi e]^ , £21 ^nd £3 . Furthermore, the EO material 




Figure 1: Relationship of laboratory coordinates {x,y, z) to the conventional electro-optic coordi- 
nates (1, 2, 3). SWP is taken to occur parallel to the x axis with phase velocity v_. 

may belong to one of 20 crystallographic classes of point group symmetry, in accordance with the 
relative values of the EO coefficients rjK ■ 
The rotation matrix 

(cos-0 — sin-0 \ 
sinV' cosi; (4) 

1 J 

in Equation (1) denotes a rotation about the z axis by an angle tp G [0, 27r). The matrix 



Ryix) = 



-sinx 





-1 


cosx 



cosx 





sinx 



(5) 



involves the angle x G [0, 7r/2] with respect to the x axis in the xz plane, and combines a rotation 
as well as an inversion. The angles ip and x delineate the orientation of the EO material in the 
laboratory coordinate system, the full transformation from laboratory coordinates {x, y, z) to those 
used conventionally for EO materials (1,2,3) being illustrated in Figure 1. 

Now epE is a symmetric matrix, regardless of the magnitude and direction of _E ^. With the 
assumption that all of its elements are real- valued, epE can be written as ai {UxU^ + UyUy + Uzliz)~^ 
Q^2 {UmU.n ~^ liinliim) 1 where ai and a2 are scalars and the unit vectors u^ and u„ are parallel to the 
crystallographic axes or the optical ray axes of the EO material in the (1, 2, 3) coordinate system 
(Chen 1983). Likewise, ep^ can be written as /3i {u,j.u^ + UyHy +1kzlLz) + /32 (%«„ +u„w„), where 
j3i and /32 are scalars and the unit vectors Up and u are parallel to the optic axes of the EO 
material in the (1,2,3) coordinate system. Thus, the application of the uniform dc electric field 
not only changes the eigenvalues of ep^;, but also rotates both optic axes and both optical ray 
axes, in general. 



2.1 Field representations 

Without loss of generality we assume that the SWP direction is parallel to the x axis. The fields 
in the half-space z < must satisfy the equations 






2. . ' (6) 



where Eq and ^o are the permittivity and permeability of free space, and £_g and Tig are the 
complex-valued amplitudes of the electric and magnetic field phasors in the isotropic dielectric 
material. The wave vector 

ks^ko{>iu^-iqsU^), (7) 

where x and 

qs = +y^?^^ (8) 

are the normalized propagation constant and the decay constant, respectively, whereas ko = 
bj^eoHo is the free-space wavenumber. We must have Re [qs] > for SWP; furthermore, x must 
be real-valued and positive for un-attenuated propagation along the x axis. Accordingly, the 
acceptable solutions of Equations (6) are 

S.S = ^si Uy + ^s2 {iqs Ux + ^i^z) (9) 

and 

tLs^ \ — [^sl {iqs Ux + XUz)^ ^s2 njuy] , (10) 

V Mo 

where Asi and As2 are unknown scalars. If qs is purely real-valued, the electric and magnetic field 
phasors of the surface wave decay exponentially with respect to z as z ^ — oo; otherwise, the decay 
is damped sinusoidal. 

The fields in the half-space z > must be solutions of the equations 

fee X^e =^Ato2ie 1 ,^^. 

The wave vector 

fcg = ko {xu^ + iqe Uz) (12) 

must have Re [qe] > for localization of energy to the bimaterial interface. A purely real-valued 
<7e indicates an exponential decay of the field phasors with respect to z as z ^ +00, whereas 
Im [qe] y^ indicates a damped-sinusoidal decay. 

Substitution of feg into Equations (11) yields a set of six homogenous equations that arc linear 
in the six Cartesian components of £g and Tie- Setting the determinant of the coefficients equal to 
zero gives the dispersion equation 

Co + Ciqe + C2ql + C^ql + dqt = , (13) 



where the coefficients 



2 1^ 2 !^ 2 c\ 

— ^xx^yz ^ ^xy^^^ ' ^xz^yy ^xx^yy^zz ^^xy^xz^yz 

./ I _2_2\2_ 4 

~rl,^2:a;^yy ~r *^2:a;^j22: ^xy ^xz)^ ^xx^ i 

Ci = ^ly^xz^yy ^ ^xy^yz)^ ^ ^l^xz^ j 

C2 = £3:2 + <^yz - {<^xx + <^yy)<^zz + (e^a; + ^zz)^ , 



c. 



2iKe.j 



Ci = -ez 



(14) 
(15) 
(16) 
(17) 
(18) 



involve e^y = u^ ' ^rei ■ Uy, etc. 

The solution of Equation (13) leads to four values of qe- We select the two values of qe that 
conform to the restriction Re [qe] > 0, and label them as qei and qe2- The corresponding wave 
vectors in the EO material are denoted by k^^ and k^2- Accordingly, with unknown coefficients 
Aei and Ae2^ we write the Cartesian components of the field phasors in the half-space z > as 



Ifec ' tie 



Uy -Se 



ILz ' Slc 



^xyf-yz + £a;zgel + jj^yylel + gel)^ ^ ^i^xz + gel)>^^ 



Ael 



- ^xyf^yz + ea;z9e2 + i{<^yy<le2 + 9e2)^ ~ ii<^xz + (le2)>t^ 



A, 



'^xy'^xz I '^xx^yz 



<^yz<le2 - ^exyqe2^ 



Ael + Ae2 , 

{exxeyy - 4y + {e^^ + eyy)q^^ + qj-^) - (e^^ + gg^);^^ 

^xx^yz ~r ^xy^xz ^yzQel ' ^^xyQel^ 

(2 /' \2 4^/' 2\2 

^xx^yy — e^y + (ea;a: + ^yy)Qe2 + ?e2J ^ (^2:a^ + ?e2)^ 

^xx^yz ~r ^xy^xz ^yzQe2 ~^ ^^xyQe2^ 



Ael 



A 



e2 



and 



ILx -He 



Uy-lLe 



= I 



-i. —qel{Ael + Ael), 
V Mo 

'^yz I ^xz^yy jQel t ^xzQei ~r (^ "^a;a; r ^xx^yy ^^^y 




«(- 



^xy^yz 



eljK 



^xx^yz ^xy^xz ~r ^yzQd ^^xyQel^ 



\ ^xx I ^xx^yy ^xy ^^xzQel ^^ ^xxQel J-^ 



■'-e.xx>^^ 



+H — 
Mo 



^xx^yz ^xy^xz H~ ^yzQ^i ^^xyQel 

^ ^{-exyeyz + exzeyy)qe2+exzql2 



■A 



-el 



+ ( ^xx + ^xx^yy ^xy)-'^ 



^xx^yz ^xy^xz i ^yzQe2 ^^xyQe2^ 



\ ^xx i ^xx^yy ^xy ^^xzQe2 ' ^xxQ_e2}^ ^xx^ I . 
9 '■ ) Ae2 



^xx^yz ^xy^xz i ^yzQe2 ^^xyQe2^ 



-e2 , 



Uz-lLe 



' — >c{Ael+Ae2). 

^J'o 



(19) 
(20) 



(21) 

(22) 
(23) 
(24) 

(25) 

(26) 
(27) 



2.2 Boundary conditions 

The boundary conditions at the interface z = lead to the following four equations: 

■ Ti.„ — 



Uy ■ £, 



Uy-Us = 




(28) 



These equations may be cast in matrix form as 



M ■ 



( ^^1 \ 

Ael 
\Ae2 ) 








(29) 



where M is a 4x4 matrix. For a non-trivial solution, the determinant of M must equal zero; thus, 
the SWP dispersion equation is 

dot M = . (30) 

Because of the complexity of Equation (30), an algebraic result could not be obtained and recourse 



was taken to a numerical method of solution. Parenthetically, when \E_ 



0, the Pockels effect 



is not invoked, and the presented formulation simplifies to that of Polo et al. (2007a) 



3 Numerical Results and Discussion 



The Pockels effect occurs only in dielectric materials that lack inversion symmetry. Some examples 
of materials of this type are potassium niobate, lithium niobate and gallium arsenide. We chose 
potassium niobate for calculations since it has very large EO coefficients. As potassium niobate 
belongs to the orthorhombic mm2 class, the only non-zero EO coefficients are ri3, r23, r^s, r42, 
and rsi; hence, we get 



/ .r(i 



£PE 



V 



1 


(0)J0) 



e[°^ nsEt) 







4°^(i-4°^23i?r 



r^iEf 



'dc 



^2 



_JO)JO) 



\ 



7-42 -B; 



r42E: 



]dc 



.(0) 



(1 



.(0) 



(31) 



^33^3^^) / 



from Equation (2). Constitutive data for potassium niobate are as follows (Zgonik et al., 1993): 



.(0) 



^33 



ci - 4.72, e^"^ = 5.20, e^'^^ = 5.43, rig = 34 x lO^^^ m V"!, rag = 6 x lO^^^ m V"!, 
63.4 X 10^12 j^ V"\ r42 == 450 x lO^^^ ^^ V"\ and r^i = 120 x lO^^^ ^^ y^^ 

The existence of a surface wave was determined at various values of x, "07 s-nd n^ by satisfactory 
solution of Equation (28). Purely for illustrative purposes, the most detailed calculations were 
performed at only one value of x- 60°; at this value of x the range of propagation directions is 
relatively wide. Let us note that the range of values of Ug for SWP at the planar interface of a 
non-EO biaxial dielectric material and a non-EO isotropic material is limited (Polo et al., 2007), 
and the range of values of the orientation angle ij) for a specific Us is quite small as well. In the 
remainder of this section, the mid-point of the 7/;-range is denoted by '0™ and the width of that 
range by A^. 
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Figure 2: Median propagation angle ipm and the width A-ip of the V'-range, when x 
\E \ = 0. See the text for other parameters. 



60° and 



3.1 Results for x = 60° 

In order to provide a basehne for the influence of the Pockels effect on SWP, Figure 2 shows tpm 
and Aijj as functions of rig for x = 60° when \E '^\ = 0. The plots in the figure cover the entire 
ris-range over which SWP was found possible. The n^-range extends from approximately 2.292 to 
2.33, and is thus only 0.038 in width. The ^-rangc in the figure is limited to [0°, 180°]; if SWP is 
possible for a certain value of tp, it is also possible for —ip, when \E'^'^\ — 0. 

Over the range 0° < -0 < 180° two bands of tp values for SWP can be deduced from Figure 2 
for each value of Us, except at the largest value of n^ where the two bands coalesce. At the lower 
limit of the rig-range, ipm approaches either 0° or 180°; while at the upper limit of that range, tpm 
approaches 90° for both bands. When one considers the entire ^-range (i.e., —180° < ip < 180°), 
there are four bands of ^'-values that merge into two bands at both limits of the rig-range. 

In Figure 2b, Aip is shown as a function of rig. Only one curve is shown since both r/^-bands have 
the same width at each value of Ug. The curve has a single peak and A-0 goes to zero at the two 
endpoints of the rig-range. The maximum value of A^' is about 0.03° and occurs at rig w 2.296, 
i.e., near the lower endpoint of the rig-range. Thus, the ^-bands for SWP are really narrow. 

In order to explore the influence of the Pockels effect on SWP, \Ef'^\ was next set equal to 
10^ V m^^. This dc field was oriented along each of the laboratory coordinate axes {x,y, and z) 
separately. We now describe the results for each direction of orientation of the dc electric field in 
order of increasing complexity. 

Let us begin with E ^ ^ u^ 10^ V m^^. Both ipm and Aip are plotted against rig in Figure 3. 
Just as for the plots for |£^ '^j = already discussed, there are four f/^-bands with mirror symmetry 
about the rig-axis. Figure 3a shows that the rig-range for SWP has grown slightly and shifted to 
lower values of rig compared to Figure 2a. With approximate lower and upper endpoints of 2.286 
and 2.327 respectively, the width of the rig-range is 0.041. Both i/j-bands meet at V' = 90°, which 
is just the same as in the absence of the dc electric field. In Figure 3b, a single curve describes 
the widths of both V'-bands as rig varies, which is similar to the curve in Figure 2b for \E_ '^\ =0. 
The height of the Atp-peak, however, is approximately 0.062°, a little more than double that of 
the peak for \E'^'^\ ~ 0, and occurs at Ug « 2.289. 

Figure 4 shows the ipm-ns and A-0-ng curves when E_ '^ = u^ 10^ V m^^. Just as in the previous 
two cases, there are four i/i-bands with mirror symmetry about the rig-axis. So, only two i/'-bands 
are shown in the figure. Evidently from Figure 4a, the rig-range for SWP is larger than when the 
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Figure 4: Same as Figure 2, except that E_ ^ = ILx 10^ V m 



dc electric field is either absent or aligned parallel to the z axis. The width of the ris-range for 
SWP has grown to 0.049, with the lower and upper endpoints of the rig-range being 2.289 and 
2.338, respectively. Although both V'-bands in Figure 4a have wider ng-ranges than in the previous 
two figures, the upper band (labeled Band 2), has increased more than the lower band (Band 1). 
In addition, the two t/j-hands meet at a lower value of ipm (~ 60°), and, thus, lack the symmetry 
about tpm — 90° seen in Figures 2a and 3a for the cases oi E_ '^ = and E_'^ \\ Uzi respectively. 

Both ip-hands in Figure 4b still show a single Af/j-peak each, towards the lower endpoint of the 
rig-range. However, the A-ip-ris curves for the two bands are not identical. Whereas the A-^-peak 
for Band 2 is about the same as for E '^ \\ u^ in Figure 3b, the Aip-peak for Band 1 is about a third 
of that value. In addition, the positions of the Atp-peaks are shifted slightly from the value for 
\E_'^''\ = 0, with the peak for Band 2 shifted downward to rig ~ 2.300 and that for Band 1 upward 
to Us = 2.292. 

Figure 5 displays the influence of the dc electric field when applied along the y axis, i.e., E_ "^ = 
u 10^ V m^^. In Figure 5a, the full range of tp, [—180°, 180°], is displayed as the i/;-bands for SWP 
are no longer symmetric about V' = 0°; instead, these bands (labeled I to IV) are now located 
symmetrically about r/; = 90°. The widths A^ of the ^-bands are shown in Figure 5b. As is the 
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Figure 5: Same as Figure 2, except that E_ '^ ~ Hy 10^ V 



case for the dc electric field applied parallel to the x axis, the widths of all four bands are not equal. 
Band I with a i/'-range of approximately [—27°, 17°] and Band IV with an approximate ^-range of 
[163°, 233°] have the same width as a function of n,,. The maximum width of these two i/'-bands 
is about 0.052° and occurs near Ug — 2.304. Similarly, Band II with a ■0-range of [20°, 84°] and 
Band III with a V'-range of [96°, 160°] share a common Atp-Ug curve. The widths of Bands II and 
III are more than a factor of 10 smaller than of Bands I and IV, the A-0-peak for Bands II and III 
being 0.0028° at n^ « 2.299. 

The sudden change in At/j at the upper endpoint of the Us range for Bands I and IV in Figure 5b 
should be noted. A similar sudden change was also found for x — 75°. The region of coalescence 
of two bands on the rig-axis is hard to delineate accurately. Quite possibly, the sudden change is 
a numerical artifact; it will be the subject of future investigations. 

The effect of variation of the magnitude of the dc electric field is shown in Figure 6 with a plot 
of tpm vs. the dc field's signed magnitude, for E_ '^ parallel to the x, y, and z axes. For the dc field 
oriented along both the y and the z axes, V'm increases as the signed magnitude becomes more 
positive: the relationship is nearly linear for _E "^ || u^, and somewhat S-shaped for ^ '^ \\ u . On 

the other hand, when ^ "^ |j u^., f/'m decreases as the signed magnitude becomes more positive. The 
foregoing results clearly show that by varying the magnitude and/or direction of the dc electric 
field, the direction of SWF, relative to the crystal axes of the EO material, can be controlled. 



3.2 Results for other x 

Limited results were also obtained for x = 0° , 30° , and 75° , with the search for SWP restricted to 
V- e [0°,90°]. Figure 7 contains Vm-n^ curves for x e {0°, 30°, 60°, 75°}, for (a) Ef" = and (b) 
Ef''^ of signed magnitude +10^ V m~^ and parallel to the x, y, and z axes. Although both Band 
1 and Band 2 (when E'^'^ \\ ILx) and Band I and Band II (when E'^'^ \\ Uy) should be visible in the 
plots, for simplicity, only Band 1 and Band II are shown. 

The curves for the chosen values of x ^^re isomorphic. At each value of x, a marked downward 
shift of the ipm-ris curve, compared to the case for \E_ '^\ — 0, occurs when ^ "^ || u^., and an 
upward shift for E'^" \\ - ^^ — - "'— " "^'^'^ ^"" ^"^^ 



There is also a shift for E' 



but it is much smaller and is almost 



unnoticeable at the lower values of x; however, the shift is significant for x — 75°, on the same 
order as when the dc electric field is applied along the other two axes. 

Figure 8 shows Atp-Us curves for x '= {0°, 30°, 60°, 75°}, for the dc electric field configured as 
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Figure 6: Median propagation angle tpm vs. the signed magnitude of _E ^ for E_ '^ parallel to the 
X, y, and z axes; Us — 2.32, x — 60°. See the text for other parameters. 

for Figure 7. When x = 0°, the curves for all three orientations of the dc field and the zero dc field 
are similar, having nearly the same peak values of A^ and peaking at close to the same values of 
Us- As X increases, the curves become differentiated both in the Ai/i peak height and its position 
on the Us axis. Among the four values of x explored, the largest peak values of A^ occur for 
X = 60°. This is particularly true when E_ '^ \\ u^; then some values of Atp are more than an order 
of magnitude larger than found at other values of x- 



4 Concluding remarks 

The influence of the Pockels effect on SWP at the interface between potassium niobate and an 
isotropic dielectric material has been demonstrated by our numerical studies. With the application 
of a dc electric field, we have noted the shift in the n^-range and the V^-bands that permit SWP. 
The greatest influence of the Pockels effect is on the median propagation angle tpm of the very 
narrow f/^-bands in which SWP is possible. Shifts of over 30° at a fixed value of rig have been 
deduced. Thus, the Pockels effect can be pressed into service for an electrically controlled on-off 
switch for surface waves. 

As of now, the existence of surface waves at the planar interface between a biaxial dielectric 
material and an isotropic dieletric material has not been demonstrated experimentally. The narrow 
■0- and ris-regimes for SWP may discourage experimentalists from searching for surface waves in 
these scenarios. Electrical control of SWP direction may provide a convenient way to search for 
surface waves. The experiment could be carried out with a fixed geometry defining the propagation 
direction and the orientation of the linear EO material, and the signed magnitude of E_ '^ could 
then easily be swept electronically until a surface wave is detected. Various experimental configu- 
rations for exciting surface waves are available (Agranovich and Mills, 1982), as also are optically 
transparent electrodes to apply dc electric fields (Minami 2005; Medvedeva 2007). Finally, we hope 
that our work shall spur the development of artificial linear EO materials with much higher EO 
coefficients than now available. 
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Figure 7: Median propagation angle i/'m versus Ug for: a) x == 0°, b) x = 30°, c) x = 60°, d) 



X = 75°. Only Band 1 is shown for S '^ \\ u^; only Band II is shown for E_ '^ \\ Uy. See the text for 
other parameters. 
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